Triplet spin resonance of the Haldane compound with interchain couphng 
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Spin resonance absorption of the triplet excitations is studied experimentally in the Haldane mag- 
net PbNi2V208. The spectrum has features of spin 5 = 1 resonance in a crystal field, with all three 
components, corresponding to transitions | ± 1) t:> |0) and | — 1) ^ being observable. The 
resonance field is temperature dependent, indicating the renormalization of excitation spectrum in 
interaction between the triplets. Magnetic resonance frequencies and critical fields of the magnetiza- 
tion curve are consistent with a boson version of the macroscopic field theoryj^ii^ implying the field 
induced ordering at the critical field, while contradict the previously used approach of noninteracting 
spin chains. 
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I. INTRODUCTION 



The inorganic dielectric PbNi2V208, with the mag- 
netic structure formed by chains of Ni^+ ('5'=1) ions, ex- 
hibits a Haldane like ground st&teJ^ Due to anisotropy 
and interchain exchange the Haldane energy gap is re- 
duced and PbNi2V208 is close to the critical point of 
the quantum phase transition from a spin-liquid to an 
ordered easy-axis antiferromagneti^ Besides, the spin- 
liquid phase may become unstable at the critical mag- 
netic field, corresponding to the vanishing of the energy 
gap of triplet excitations. In Heisenberg exchange ap- 
proximation the critical field value He and the energy gap 
A are related by a simple relation gfisfJ-oHc = A. The 
influence of the single-ion anisotropy on the Haldane spin 
chains was analysed theoretically by the exact diagonal- 
ization for finite chains^ by a perturbative approach)^ 
as well as by macroscopic field theory methods<^ In the 
magnetic field range far below He all models result in the 
same energy levels, parametrized by the main gap, one or 
two anisotropy constants and (/-factor. The perturbative 
description of excitations in the Haldane chain appears 
to be identical to the description of an isolated spin 5 = 1 
in a crystal fieldJ. Extrapolated to the field of spin-gap 
closing, the perturbative approach yields critical fields 



gUBlJ'oHac = v^A^A^, 



(1) 



here the magnetic field is applied along principal direc- 
tion a, while (3 and 7 note other principal directions. For 
the uniaxial case the Sz — ±1 triplet components with 
the momentum of tt have equal gaps = Aj^, and the 
gap of 5*2 = component, Az, is higher or lower depend- 
ing on the anisotropy type. The value of the gap splitting 
caused by the crystal field is, in terms of 5 = 1 problem, 
the effective anisotropy constant -De// = A^ — A^. For 
Haldane chains De// = -1.98i:> and {2Ax + Az)/3) = A, 



where D is a single ion anisotropy constant contributing 
to the single ion Hamiltonian with the term DS'^. 

The spectrum of a real spin-gap system may be 
more complicated, as, e.g., for the Haldane magnet 
Ni(C2H8N2)2N02(C104) (abbreviated as NENP), when 
the spin gap remains non-zero due to the staggered g- 
tensor;^ or for the dimer spin-gap system TlCuClg^iiS 
with a nonlinear frequency-field dependence near the 
critical field, at which the field-induced 3D ordering oc- 
curs, etc. PbNi2V208 is a convenient model system for 
study the influence of 3D correlations and anisotropy on 
the spin-gap in a magnetic field, because this compound 
demonstrates a spin-gap split by the anisotropy, spin- 
liquid behavior below the critical field and field-induced 
3D antiferromagnetic ordering above H^i^ 

The perturbative approach of noninteracting Hal- 
dane chains (in particular, formula ([1])) was usedi^ for 
PbNi2V208 to derive the energy gaps from the crit- 
ical field values. Further, the velocity of spin excita- 
tions, intrachain and interchain exchange integrals were 
found by fitting the inelastic neutron scattering inten- 
sity from a powder sample keeping the above spin gap 
values fixediii^ii^ Nevertheless, the extrapolation of the 
perturbative approach to high fields may be doubtful. 
While it is justified by the exact diagonalization of the 
finite chain problem;^ the field-theory treatment^ has re- 
sulted in two models, diverging in high fields. The boson 
model should be exact in presence of interchain coupling 
of nearly a critical value, while fermion model is exact 
for a hypothetic Hamiltonian. Analogous theory consid- 
ering a general case of a spin-gap magnet with the field- 
induced 3D ordering was developed in Ref The source 
Lagrangian and relations between the critical fields in the 
bosonic model^ and in the macroscopic theorjsi^ are iden- 
tical. The results of the fermionic model are close to the 
perturbative theory and exact diagonalization of finite 
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FIG. 1: Pulse field magnetization of the aligned sample of 
PbNi2V208. Solid curves; experiment, dashed stright lines - 
extrapolations to j^cx,| 

chain problem, while the boson model and the modeliS 
diverge from the perturbative curves near the critical 
field. Probably, this divergency should be ascribed to 
fluctuations, which are suggested to be suppressed (e.g. 
by 3D coupling) in the boson model of Ref. m and in the 
macroscopic theory of Ref. ITol . 

In the present paper we study the triplet excitations 
by means of the high-frequency electron spin resonance 
(ESR) spectroscopy which enables one to measure di- 
rectly the frequencies of the transitions between the spin 
sublcvels of the collective triplet states. We recover the 
splitting of the zero-fleld energy gap from the experi- 
ment and relate it to the difference between the criti- 
cal fields Hxc and H^c- The zero- field splitting of the 
triplet levels appears to be about one-half of the value 
estimated by the perturbative relation ([T]). At the same 
time, the spectrum of triplet excitations is well described 
by macroscopic models, implying 3D ordering at the criti- 
cal fieldi^ii^ Thus, the influence of the interchain coupling 
was found to be of importance for the magnetic excita- 
tion spectrum in the vicinity of the critical field. 

II. EXPERIMENT 

We used ceramic samples of PbNiV208 prepared as 
described in Reffisl and magnetic field aligned samples 
prepared following Ref. [l| with the orientation of c-axis 
of crystallites parallel to the aligning field. Magnetiza- 
tion curves taken by the pulse technique are shown in 
Fig. [T] The critical field values derived from these curves 
are /ioi?c||=19 ±0.5 T and /ioffc-L=14±0.5 T which are 
well consistent with the data of Refslll and[ll|. The ESR 
spectra were measured both for ceramic and aligned sam- 
ples, as microwave transmission vs magnetic field depen- 
dences. Magnetic field was created by a 16 T cryomagnet. 
Microwave signal in the frequency range 70-500 GHz was 
generated and detected by means of the vector network 
analyzer from the ABmm system at the KYOKUGEN 
Center. 

A typical ESR spectrum of a ceramic sample and its 
evolution with temperature is presented for the frequency 
/ =150 GHz on Fig. [21 Analogous spectra for the 
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FIG. 2: 150 GHz ESR spectra of a PbNi2V208 ceramic sam- 
ple measured at different temperatures. Inset: Temperature 
dependence of the ESR fields A, B, C. Lines are guide to the 
eye. 



aligned sample and the frequency of 275 GHz are pre- 
sented on Fig. [3] At the temperatures above 30 K there 
is a single ESR line, corresponding to a typical exchange- 
narrowed Ni-^+ ESR signal with ^-factor value of 2.23. At 
cooling the sample we observe the diminishing of the in- 
tensity of this signal and appearing of three ESR lines 
apart from the high temperature ESR field. These sig- 
nals are marked by letters A, B and C on Figs. [2]and[3l 
They disappear again at low temperatures, demonstrat- 
ing a temperature dependent shift of the resonance field, 
as shown on the Insets of Figures. At low temperatures 
an upturn of the resonance absorption in g-factor range 
1.9-2.2 with the maximum at at g = 2.2 occurs. This 
low-temperature signal is due to defects and impurities. 
The intensity of marked signals, increasing with temper- 
ature, indicates a thermally activated nature of this kind 
of resonance. The separated lines of thermally activated 
absorption are distinctly seen in the temperature range 
15-20 K, at lower temperatures the intensity is weak be- 
cause of the spin gap, and at higher temperatures the 
lines merge into a single line. No absorption signal cor- 
responding to the singlet-triplet transitions is observed. 

Comparison of the ESR absorption of the ceramic and 
aligned samples at close frequencies is given on Fig. U) 
Dashed arrows near the curve, corresponding to the ce- 
ramic sample, mark boundaries of the absorption band 
and letters A, B, C denote thermally activated lines for 
both samples. The central peak is ascribed to defects 
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FIG. 3: 275 GHz ESR spectra of a PbNi2V208 aligned sam- 
ple measured at diflterent temperatures for H perpendicular 
to the alignment axis. Inset; Temperature dependence of the 
ESR fields A, B, C. Lines are guide to the eye. 
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FIG. 5: Frequency-field diagram of the thermally activated 
magnetic resonance of PbNi2V20g ceramic and aligned sam- 
ples. Triangles and squares present experimental data for 
resonances marked by A, B, C. Lines are model calculations. 
Solid circles represent the extrapolations to T=G noted by 
curved arrows. 



and diminishes with temperature as a conventional para- 
magnetic resonance line. The analogous record of the 
aligned sample at i? _L c has more distinct resonance 
lines at the positions B and C, than the ceramic sample. 
The positions of the local maxima of thermally activated 
absorption, taken at different frequencies, are plotted on 
Fig. El 



III. DISCUSSION 
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FIG. 4: 300 GHz ESR spectrum of a PbNi2V208 ceramic 
sample at T =20 K (lower curve) and 275 GHZ spectrum of 
an aligned sample (upper curve). 



A. General aspects 

The lineshape of the thermally activated absorption 
(Figs. [2]and[l]) corresponds approximately to the known 
ESR spectrum of a powder and crystal samples with 
5 = 1 ions in an axial crystal field (see, e.g., Ref. 0). 
For a uniaxial crystal, this ESR spectrum, in case of 
gfJ-BfJ-oH ^ D, consists of three resonance lines corre- 
sponding to the transitions | ±1) ^| 0) with ASz = ±1 
and I —1) ^1 1) with ASz = ±2. Two resonance lines, 
corresponding to ASz = il, have larger intensity and 
are on both sides from the free spins resonance field 
^o^^o = 27r?i//(g/LtB). The distance between these two 
lines is 2D/giJB for H \\ c and approximately D/gfiB 
for H J- c. Here D is a single ion anisotropy constant 
contributing to the single ion Hamiltonian with the term 
DS^. Naturally, for a powder sample there should be a 
band of absorption instead of resonance lines, the whole 
band width is 2D/g^B, and the maxima of absorption 
are near the resonance field of crystallites with c J- H 
because of their largest statistical weight (the field inter- 
val between maxima is again D/gfis)- In case of the easy 
plane anisotropy the intensity of the lower field maximum 
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FIG. 6: Triplet energy levels in a spin gap magnet with 
fioHc±=14:.0 T and /^o^^c|| = 19.0 T, calculated for noninter- 
acting Haldane chains^ (thin lines) and derived in the macro- 
scopic modeli^ (thick lines). Vertical segments present ob- 
served transitions at /=300 GHz 



is larger than that at the upper one because of a higher 
population of the lower spin sublevels. For the case of the 
easy axis anisotropy the relation between the intensities 
should be inversed. 

Applying the above consideration of 5 = 1 magnetic 
ion in a crystal field we conclude, that the effective 
anisotropy of the triplet excitations is of the easy plane 
type {Deff > 0). According to Ref. 0, -De// and D are 
of different signes, thus, the single-ion constant D should 
be negative, which agrees with the easy axis anisotropy 
of the impurity induced ordered phasei^^ 

Maximums of the absorption are separated by approx- 
imately 3.5 T, which corresponds to D^ff — 110 GHz 
(here and farther on we assume g = 2.23, as proved by 
ESR at T > 30K). The leftmost absorption line demon- 
strates a characteristic frequency-field dependence (Fig. 
[5]) with asymptote / = 2gfiBfJ-oH/{2TTh), marking a two- 
quantum transition with ASz = ±2. Its zero-field gap 
is equal to the splitting of the triplet sublevels^. Ob- 
servation of this absorption at the frequency 150 GHz 
indicates directly, that the splitting of the triplet levels 
by the effective crystal field is less than or about 150GHz. 

The perturbative approach^ yields by ([T]) the zero-field 
gaps E,mni = 320 ± 10 GHz and £:„„2 = 590 ± 10 GHz. 
The corresponding splitting of 270 GHz is twice as large 
as estimated from the observed ESR lines. Alternatively, 
energy levels of the triplet excitations can be calculated 
within the boson model^^ii^ which imply 3D ordering 
above the critical field. Details of the calculation are 
given in Appendix. For D^ff > the gaps are related to 
the critical fields as follows: 




(2) 



FIG. 7: Temperature dependence of the resonance field of the 
mode B at the frequency 475 GHz for the aligned sample at 
H ± c. Dashed lines are linear and second order polynom 
extrapolations to zero temperature. 



From the critical field values we get here Emini — 
440 ± 10 GHz and Ermn2 = 590 ± 10 GHz. The cor- 
responding splitting of 150 GHz is in much better agree- 
ment with the above experimental results. The calcu- 
lated triplet sublevels for perturbative and macroscopic 
models are presented at the FiglHl Here all theories are 
locked to the observed values of critical fields. The dif- 
ference in Eraini occurs due to the nonlinear field depen- 
dence of the energy of the lowest triplet sublevel near 
Hci^ in the macroscopic models. It should be noted, that 
in low fields, i.e. far from the critical point, all mod- 
els would demonstrate identical results when they would 
locked on the same energy gaps at if = 0. 

The frequency-field dependence for the case _L c 
calculated in a perturbative approach and in the macro- 
scopic models is presented in Fig. [5j It should be noted 
that no fitting parameters are used here, the value of zero 
field splitting, -De// = Emin2 — Emini is taken from crit- 
ical fields using relations (21). The resonance fields were 
taken at T =15 and 20 K, while critical fields are mea- 
sured at 4.2 K. Thus, we have to extrapolate the observed 
resonance fields to zero temperature. The temperature 
evolution of the ESR spectra, analogous to that shown 
in Fig. [31 was followed for the frequencies 150, 275 and 
475 GHz. The resonance fields, extrapolated to T=0 are 
shown on FigO by solid circles, pointed by curved ar- 
rows. The data, extrapolated to zero temperature, lie on 
the model curve within the experimental error. 

Note that the energy levels predicted in Refs. [H, d and 
[121 by use of the perturbative model^ {^eff = 270GHz 
and Emin2 / {^Trh) =590 GHz) does not correspond the 
above experimental data to any extent. Thus, the ex- 
perimental ESR frequencies are in a agreement with the 
energy levels predicted by the macroscopic theories j^^i^, 
rather than with the theory considering isolated Haldane 
chains 



B. Modeling 

We performed a model calculation of the ESR absorp- 
tion for the powder sample of a spin gap magnet with 
gapped triplet excitations (see Figs. |niand|51). The mod- 
eling was performed by two methods, giving practically 
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FIG. 8: 275 GHz ESR spectrum of a PbNi2V208 aligned 
sample at T =15 K, i/ _L c, and the calculated spectrum in 
S — 1 model for the triplets at the bottom of the excitations 
zone (q=0) and a spectrum, averaged over the excitations, 
excited at T=2Q K. Parameters used =590 GHz, 

D^ff{0)= 150 GHz, ESR linewidth 10 GHz, s=4895 GHz, 
g = 2.23 
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FIG. 9: 300 GHz ESR spectrum of a PbNi2V208 ceramic 
sample at T =20 K (black solid line), the same spectrum 
with the defects Lorentzian line subtracted and model spectra 
as described in the text. Vertical lines are at the resonance 
fields for crystallites with the field lying in the basal plane. 
Parameters used are: i?min2=590 GHz, _De//(0)= 150 GHz, 
ESR linewidth 10 GHz, s=4895 GHz, g = 2.23 



similar results for the lineshape and positions of the max- 
ima of ESR absorption. 

For the first method we considered the energy levels of 
effective spin S — 1 ascribed to a triplet excitation mode 
with a wavevector fe. The population numbers for the 
= ±1 and Sz = spin sublevcls of a triplet excitation 



mode were calculated according to the Plank's distribu- 
tion function. The effective anisotropy constant, equal 
to the zero field splitting of a triplet sublevels is defined 
as 



(3) 



here uj^^yiO, k) and lUz{0, k) are zero- field cyclic frequen- 
cies of the triplet excitations carrying effective spin pro- 
jection S"^-^^ equal to ±1 or zero correspondingly. We 
calculated the microwave absorption for effective spins 
5 = 1 in a uniaxial crystal field according the standard 
ESR procedure, i.e. calculating transition frequencies, 
matrix elements of the transitions and difference of the 
population numbers^. The spectrum of the triplet exci- 
tations in zero field is taken in the forr" ^'^°'^^ - 



(4) 
(5) 



with q = TT - ka, s=4895 GHz (Refi) , and 
w^(0)/(27r)= 440 GHz, u;^,y(0)/(27r)=590 GHz derived 
from the fields and Hc± according to Ref. [l^ as 
shown on Fig. [6l The field dependences of huja{k) are 
derived as the energies of corresponding sublevels of an 
effective spin 5 = 1 in a crystal field with the anisotropy 
constant Deff{k). 

This modeling also provides the low-field AS* = ±2 
transition, which is allowed in the canted orientations of 
the magnetic field with respect to anisotropy axis. 

For the second method we used the excitations frequen- 
cies, given by the macroscopic modelsj^ii^ These frequen- 
cies can be calculated for an arbitrary H and q (see Ap- 
pendix) in the low frequency approximation hiOa ^ J- 
Since this model does not provide information on ma- 
trix elements of the transitions, which are necessary to 
calculate intensity of the corresponding ESR absorption, 
we have used this model only at high fields taking the 
matrix elements for ASz = ±1 transitions to be equal 
and field independent. Population numbers of the exited 
states were taken according to the Plank's distribution 
function again. 

The model microwave absorption due to the excita- 
tions with the energy at the bottom of the zone (q=0) is 
shown in Fig. [8] along with the model absorption curve 
at T=15 K, when the excitations away of the zone bot- 
tom are also excited. Despite of the smaller value of 
Deff of the high energy excitations, the peaks of absorp- 
tion remain at the same positions due to a large spectral 
density at the zone bottom. At heating, the absorption 
curve becomes asymmetric with enlarging of the intensity 
in the field range, corresponding to resonance for smaller 
Deff, as predicted in Ref. 0. For modeling of the ESR 
absorption of the ceramic samples, we made the aver- 
aging over the random orientations of crystallites and, 
again, over the one-dimensional k-space. The calculated 
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lincshapc, presented by the dashed curve on Fig. [9]corre- 
sponds quahtatively to the observed hne, if we subtract 
the central peak of absorption, produced by the defects, 
as shawn on the Fig. It should be noted that the sin- 
gularity in the form of the step of absorption is close 
to the maximum of absorption for the single-fc calcula- 
tion, both singularities arc due to the maximum of the 
statistic weight of the crystallites with the four fold axis 
perpendicular to the magnetic field. For the fc-space in- 
tegrated absorption the steps are at the same positions, 
while the maxima are smeared due to the contribution 
of the higher energy excitations which have smaller value 
of Deff{k). The results of the modeling are practically 
independent on the natural ESR linewidth changing in 
the range 1-10 GHz. 



C. Temperature dependent ESR 
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V. APPENDIX 

Following Rcfflo!, the Lagrangian for the vector field 
T] corresponding to S* = 1 excitation for the crystal with 
the uniaxial symmetry is 



The temperature dependence of the triplet ESR field 
observed in the present work may be. probably, at- 
tributed to an interaction between the triplets, analo- 
gous to the spin-wave frequencies in conventional ferro- 
and antiferromagnets, which is known to be renormal- 
ized with excitation of magnons (see, e.g., Refs. [l^ and 
[iTh . The temperature dependence of the triplet excita- 
tion spectrum was observed also for the dimer spin-gap 
compound TlCuCla .— An alternative reason for the tem- 
perature dependence of D^ff and the triplet excitation 
spectrum may be the change of the correlation length 
with heating. However, according to experiments^^, the 
correlation length is reduced only for 20 % at the temper- 
ature of a half of the intrinsic Haldane gap 0.41J, as in 
our case. Apparently, this may not substantially change 
the effective anisotropy constant of triplet excitations. 
The excitation of the triplets apart from the bottom of 
the excitation spectrum also may not change the posi- 
tion of the absorption maxima, despite they should have 
smaller Dej / . As noted in Ref. 0, the influence of higher 
energy excitations will result in an asymmetric form of 
the ESR lines but not in the shift of the absorption max- 



(6) 



here 7 = guB^/fT'- 

Corresponding equation of motion has the form 



rj + 2j['nH] ~ YH^r] + YH{r]H) + Ar] 

Vx 



-G.jd.d.v b\ r/y I = 0. (7) 



The frequencies at k ~ ir/a and H = are 



,{o) 



'A.^^ = VA-b 



(8) 
(9) 



The field and fc-dependences for principal directions of 
the magnetic field are: 
for H 11 z: 



IV. CONCLUSIONS 

In conclusion, the ESR spectrum of triplet excitations 
in a Haldane magnet PbNi2V208 demonstrates a tem- 
perature dependence, probably, due to the interaction 
between the excitations. ESR of triplet excitations is 
satisfactory modeled without fitting parameters, using 
the experimental values of critical fields. The relation 
between the critical fields and energy gaps corresponds 
to the boson model of Ref. and macroscopic theoryiS 
and contradicts the fermion model of Ref. and isolated 
chain calculations of Ref. 0- This indicates, supposedly, 
a strong influence of the 3D correlations on the spectrum 
of excitations in PbNi2V2 08. 



for H II 



UJQ 



(jj. 



x,y — 
UjQ = a/a 



zz I s2q2 



(10) 

(11) 



2 Tj2 I -^xx + ^zz 



2 2 I 
S q ± 



{A^x Azz) 



2j^H^ {A,,+A,,) 



A 2 cj2 2 2 



1/2 



(12) 
(13) 



here q — k — n / a. 
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For an arbitrary orientation the frequencies may be 
found by solving the secular equation 

+ (2/1^ + A,, - A^,)x' - /i2(4(A,, + q^s") ~ 
-h" + [A,, - A.,,)[l - 3cos2 e))x - 
-{A,, - A,,)h^ cos^ 0(4(A,, + qh^) -h^)=0 (14) 

here h = -jH , x = — tj^ + q^s^ + A^-x, S is the angle 
between the magnetic field and z-axis. 

Critical field in the arbitrary orientation of the applied 
field for the effective easy-plane anisotropy of triplet ex- 
citations (6 < 0) is given by the equation 



lHc{0) = 




(15) 

while for the effective easy-axis anisotropy of triplet ex- 
citations (6 > 0) He is orientation independent 

7i?c = —4°^ (16) 
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